Vortex matter in superconducting mesoscopic disks: Structure, magnetization, and 

phase transitions 
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The dense vortex matter structure and associated magnetization are calculated for type-II super- 
conducting mesoscopic disks. The magnetization exhibits generically first-order phase transitions 
as the number of vortices changes by one and presents two well-defined regimes: A non-monotonous 
evolution of the magnitude of the magnetization jumps signals the presence of a vortex glass struc- 
ture which is separated by a second-order phase transition at Hc2 from a condensed state of vortices 
(giant vortex) where the magnitude of the jumps changes monotonously. We compare our results 
with Hall magnetometry measurements by Geim et al. (Nature 390, 259 (1997)) and claim that 
the magnetization exhibits clear traces of the presence of these vortex glass states. 

PACS numbers: 74, 74.60.Ec, 74.76.-w 



Electronic device miniaturization in semiconductors 
has recently reached the ultimate limit where the num- 
ber of electrons present in the device can be tuned at will 
even down to a single electron . These systems have 
been given the name of quantum dots or artificial atoms 
since their generic properties are determined by their few- 
electron configurations much the same as in real atoms. 
Analogies can be drawn between these artificial atoms 
and type-II superconducting mesoscopic disks in perpen- 
dicular magnetic fields where the role of the electron is 
played in this case, not by the Cooper pair, but by an- 
other fundamental entity: The vortex. When the dimen- 
sions of the disk are comparable to the coherence length 
^ only few vortices can coexist in the system. In con- 
trast to the usual triangular arrangement in bulk, com- 
plex and unique vortex structures are expected to occur 
due to the competition between surface superconductiv- 
ity and vortex-vortex interaction. When the dimensions 
of the system are much smaller than ^ the very notion of 
superconductivity needs to be revised 1^. 

Transport experiments have contributed in a decisive 
way to unveil the electronic structure of artificial atoms 
Similarly, transport measurements in individual 
mesoscopic disks gave us the first experimental evidence 
of the structure of the order parameter in these systems. 
Oscillations of the critical temperature Tc as a function 
of the external magnetic field H were correctly accounted 
for by the quantization of the angular momentum L of 
the Cooper pair wavefunction or, in other words, by tran- 
sitions between giant vortex states with a different num- 
ber of flux quanta. Close to the supercondutor-normal 
phase boundary the theoretical analysis of these transi- 
tions does not present any difficulty since it simply im- 
plies solving the linearized Ginzburg-Landau equations 



Hall magnetometry g, on the other hand, is reveal- 
ing itself as a powerful tool for obtaining information of 
the order parameter through its associated magnetization 
M away from the supercondutor-normal phase boundary. 
To date, the theoretical efforts to calculate the structure 
of the order parameter and M well into the supercon- 
ducting phase have been mostly restricted to the numer- 
ical solving of the Ginzburg-Landau equations under the 
assumption of an order parameter with a well-defined 
L |p|,p^. This is justifiable for type-I superconducting 
disks. For type-II disks, however, this assumption is no 
longer valid. More precisely, it is only expected to hold 
in the Meissner state, which is associated to an L = 
order parameter, and above Hc2 where, close to the sur- 
face of the disk, the superconductivity can survive up to 
a higher critical field Hcz [^|^- This surface order pa- 
rameter is referred to as a giant vortex or macrovortex 
§,0. For Hci < H < Hc2 it has been argued § and 
shown in numerical simulations of the time-dependent 
Ginzburg-Landau equations 11 1 that the order param- 
eter can form complex structures of single-fluxoid vor- 
tices, i.e., a "budding" Abrikosov lattice. From the more 
analytic standpoint presented in this paper the appear- 
ance of these structures can only come about if the order 
parameter does not have a well-defined L. To find and 
understand the structure of these vortex states, to cal- 
culate the magnetization associated to them, and to link 
these states to the Abrikosov lattice which emerges in the 
thermodynamic limit of infinite disks are the main goals 
of this paper. We do this by expanding the order pa- 
rameter in an appropriate basis and minimizing, analyt- 
ically to a large extent, the Ginzburg-Landau functional. 
For an increasing magnetic field we find first-order phase 
transitions whenever a vortex is added to the disk. A 
non-monotonous behavior of the magnitude of the mag- 
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netization jumps signals the presence of the glassy vortex 
structures while, at Hc2, a second-order phase transition 
constitutes the signature of the condensation of the vor- 
tices into the giant vortex state. In addition we compare 
with the magnetization measured in Al disks by Geim et 
al. in Ref. |^ and claim that the magnetization jumps 
present a non-monotonous evolution which, as we stated 
above and shown below, can be associated to the ex- 
istence of single- fluxoid vortex glassy structures. This 
would indicate that the Al disks in Ref. ||] behave like 
type-II superconductors rather than type-I, possibly due 
to the expected enhancement of the effective magnetic 
penetration length in such a geometry. 

We start from the traditional Ginzburg-Landau func- 
tional for the Gibbs free energy of the superconducting 
state 
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where G„ is the Gibbs free energy of the normal state and 
[—iTW — e*A(r)/c]^/2m* is the kinetic energy operator 
for Cooper pairs of charge e* = 2e and mass m* = 2rn 
in a vector potential A(r) which is associated with the 
magnetic induction h{r). The parameters a and (3 have 
the usual meaning [ p^ . Before proceeding any further we 
must stress a not fully appreciated fact: Even for small 
values of K 1), the magnetic induction is weakly vary- 
ing in space down to fairly low fields {H sa 0.5Hc2) [ p^ . 
This observation is very important for our purposes since, 
down to -ff « 0.5iJc2, it is a very good approximation to 
consider a uniform magnetic induction [h(r) = B] and to 
expand of the order parameter in the following way: 



1 



-iLO 



L=0 



/2n 



$L(r). 



(2) 



In this expansion Cl = \CL\e^ are complex coefficients 



and 



iLe 



^^(r) are normalized nodeless functions 



that diagonalize, in the symmetric gauge, the kinetic en- 
ergy operator appearing in Eq. |l] with eigenvalues eL{B). 
We only consider disk thicknesses smaller than the co- 
herence length so that the order parameter can be taken 
constant in the direction of the field. These eigenfunc- 
tions are subject to the boundary conditions of zero cur- 
rent through the surface and the radial part ^Lif) may 
be found numerically. (The radial unit is the magnetic 
length i — a/ e*h/cB). This expansion captures both the 
simplicity of the macrovortcx (above Hc2) when only one 
Gl is expected to be different from zero and the full com- 
plexity of the order parameter (below Hc2) when several 
L's must participate. 

Direct substitution of the expansion ^ into Eq. ^ and 
subsequent numerical minimization of the resulting ex- 
pression is a daunting task bound to fail due to the large 



number of unknown variables involved. Instead, it is key 
to consider expansions in restricted sets {ii, L2, . . . , L^} 
of few N components where Li < L2 < ■ ■ ■ < L^. The 
Gibbs free energy takes the following form for each set: 
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vlft^re Gs — Gn and a are given in units of Hl^^j^-K iy 
is the volume of the disk), €l{B) is given in units of the 
lowest Landau level energy hujc/'^ (with ujc = e* B /m*c), 
R is the radius of the disk in units of ^, and B and H are 
given in units of Hc2- The terms proportional to a con- 
tain the condensation and kinetic energy of the Cooper 
pairs. All the other terms, which are proportional to 
a^, account for the "interaction" between Cooper pairs. 
There appear four types of these terms: (i) those pro- 
portional to Il{B) = J dr r reflecting the interaction 
between Cooper pairs occupying the same quantum state 
i, (ii) those proportional to IliLj (B) = J drr$|. $| . re- 
flecting the interaction between Cooper pairs occupying 
different quantum states, and (iii) the ones proportional 
to IhiLjLkiB) = J dr r $^.$1^$^^^ and proportional to 
lLiL,LkLi{B) = J drr'^L.'^Lj^Lk'^L, which, along with 
the phases (pL , are responsible for the correlation between 
vortices. The non-linear dependence on B of these inte- 
grals [as well as that of ei(i?)] comes from the existence 
of the disk surface. 

In order to find the minimum Gibbs free energy for a 
given set we have to minimize with respect to the moduli 
|CiJ, . . . , |Ci„|, the phases cpLi, ■ ■ ■ ,4'Lm of t^ie coeffi- 
cients, and with respect to B. The minimum-energy set 
of components is picked up at the end. The advantage 
of doing this selective minimization resides in our expec- 
tation that a small number of components will suffice to 
describe the order parameter for the disk sizes consid- 
ered in the experiment of Ref. ^ . As an illustrative and 
relevant example we consider in the detail the solution 
with a single component {L}. The energy functional is 
invariant with respect to the phase of the only coefficient 
so one can minimize analytically with respect to \Gl\'^ to 
obtain 
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Finally, the minimal value of B and the minimum Gibbs 
free energy for each L must be found numerically. It 
is important to notice that a disappears from the final 
expression in Eq. which leaves us with k as the only 
adjustable parameter when comparing with experiments. 
(This is also true for the more complex cases discussed 
below). The 2-component {Li, L2} solutions can be dealt 
with in a similar way. The energy functional is invariant 
with respect to the phases so one can minimize analyti- 
cally with respect to |CljP and ICl^P and numerically 
only with respect to B. The final solutions look gener- 
ically like an {L2 — Li)-vortex ring. For 3-component 
solutions (two vortex rings) the energy functional is still 
invariant with respect to all the three phases whenever 
Li + L3 7^ 2L2, and, once again, the minimization with 
respect to |CliP, ICl^P, and jCigp can be done analyt- 
ically. However, if Li+ L3 — 2L2, the two rings have the 
same number of vortices and their relative angular posi- 
tions come into play through the term depending on the 
phases. There is, however, an obvious choice for these 
phases: (fiL^ = 0, = 0, ^l, = vr. This choice gives a 
negative contribution to the free energy which reflects the 
lock-in position between the vortex rings. One important 
fact should be noted now: The components in which the 
bulk Abrikosov lattice needs to be expanded in the sym- 
metric gauge are strongly overlapping. Depending on the 
chosen symmetry of the lattice the set of components is 
either {1,7, 13, 19, . . .} for the Cg symmetry with a vortex 
at the origin or {0, 3, 6, 9, . . .} for the C3 symmetry with 
the center of a vortex triangle at the origin |l^]. Conse- 
quently, the minimum-energy solutions for disks are ex- 
pected to have strongly overlapping components which 
invalidates any perturbative treatment of the terms that 
contain the phases p^ . Moreover, unlike simpler geome- 
tries ijisj , there is no direct connection between number of 
components and number of vortices. This prompts us to 
seek solutions through numerical minimization with re- 
spect to the moduli and B for the 3-component cases just 
mentioned, and, for > 3, with respect to the moduli, 
the phases, and B whenever the terms involving phases 
are present. Fortunately, for the disk sizes like the ones 
used in the experiment of Ref. |^ we will see below that 2 
and 3-component solutions suffice to capture the relevant 
physics. 

Figure ^ shows the magnetization as a function of H 
for a disk of i? = 8^ and k = 3. (There is nothing special 
about these parameters, the only purpose of which be- 
ing to be convenient for the discussion of the physics we 
want to bring up). As indicated in the figure, different 
types of lines correspond to the magnetization obtained 
expanding the order parameter with up to A^ = 1,2,3, 
and 4 components. A common feature to all curves is 
that the magnetization exhibits many first-order transi- 
tions. Above Hc2 the A^ = 1 solutions suffice to describe 
entirely the order parameter and magnetization. Here 
the solutions correspond to a giant vortex which con- 



tains L (the quantum number of the single component) 
fluxoids. Whenever L changes by one the magnetiza- 
tion presents a (non-quantized) jump whose magnitude 
evolves monotonously with L. Below Hc2 we see that the 
N = 1 solutions underestimate the correct value of the 
magnetization. Allowing more components in the expan- 
sion has a fundamental effect: It splits the giant vortex 
into a complex structure of many single-fluxoid vortices 
[for an example see Figs. ||(a) to (d)]. This reflects in 
the magnetization curves by changing the regular evolu- 
tion of the magnitude of the jumps into an irregular one. 
Whenever a vortex is added or removed from the disk the 
symmetry of the new vortex configuration is expected to 
change which, in turn, expels the field in a different way. 
There exist configuration switches for a given number of 
vortices, but these changes do not reflect in the magneti- 
zation, in contrast to what has been suggested [^. When 
the number of vortices is large enough, there are usually 
two possible symmetries within each magnetization step: 
One with a vortex at the center of the disk and one with- 
out it which are reminiscent of the C3 and Cq symmetries 
of the regular vortex lattice. On top of the many first- 
order transitions the overall slope in the magnetization 
clearly changes at Hc2, i.e., at the transition between the 
giant vortex and the vortex glass structures. This transi- 
tion is reminiscent of the second-order transition at Hc2 
for bulk samples where M vanishes. 

Finally, we would like to point out that the magnetiza- 
tion measured by Geim et al. in Rcf. [|j presents features 
that are in good agreement with our results despite of the 
fact that the disks are made out of a strong type-I ma- 
terial as Al. In these materials the plate geometry can 
lead to vortex structures being more favorable than do- 
mains [ p6| . To compare with the experiment, we simulate 
this fact by using a higher value of k than the nominal 
one. In Fig. ^ we show the data for a disk of nominal 
radius R = 5^, thickness d = 0.6^, and k = 0.24. We 
have obtained a reasonable good agreement in the num- 
ber and magnitude of the jumps, and overall shape of the 
magnetization using i? « 5^ and k « 1 (the dotted line 
is a good example). This is consistent with an effective 
penetration length longer than expected and, possibly, 
with a coherence length shorter than the bulk nominal 
one. Although, due to the smallness of the disk, it is 
difficult to point at a second-order phase transition, the 
non-monotonous evolution of the magnitude of the mag- 
netization jumps is notorious over a large range of fields 
which, as we have shown, is an indication of the forma- 
tion of vortex glass structures. In our approximation the 
magnetic induction is uniform in space, but, given the 
good agreement with the experimental curve, this seems 
to be a much less important restriction than considering 
an order parameter with a well-defined quantum number 
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FIG. 1. Magnetization as a function of H for a disk of 
radius R — 8S, and k = 3. Different line types correspond to 
different number of components allowed in the minimization. 
For this especific set of parameters the plot shows clearly the 
necessity of considering more than one component below Hc2- 
However, no appreciable difference can be seen between the 
traces obtained using an expansion with up to A*' = 3 and 
N = 4: components. 




FIG. 2. Modulus square of the order parameter for an 
R — 8(,, K = 3 disk at H = 0.6Hc2 expanded with up to 
(a) iV = 1 (minimum-energy set {12}), (b) iV = 2 ({3, 12}), 
(c) iV = 3 ({0,3,12}), and (d) = 4 ({0,3,6,12}) compo- 
nents. The vortex structure is strongly modified from N = 1 
to N = 4 although the total number of vortices is always given 
by the largest L which does not depend on the the number 
of components considered. (Only the internal arrangement 
of vortices does.) Allowing 4-component solutions does not 
change appreciably the magnetization obtained with = 3 
(or even A = 2) at any value of H (see Fig. ^). However, the 
order parameter is still modified as can be seen by comparing 
(c) and (d). The addition of a forth component does not in- 
crease the number of vortex rings, but can help fix the relative 
position of the two existing ones in the N — 3 solution. This 
modification does not have measurable consequences, but il- 
lustrates how one of the symmetries of the triangular vortex 
lattice (the C3) can emerge, bearing confidence in our calcu- 
lations. 
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FIG. 3. Comparison between experimental data (repro- 
duced from Refs. [10]) and theory using R — 5.25^ and 
K = 1.2^. Similar considerations as in Refs. [10] have been 
followed for the adjustment of the theoretical curve. 
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